arXiv: 1503.00025v2 [cond-mat.dis-nn] 27 Nov 2015 


Metal-Insulator-Transition in a Weakly interacting Disordered Electron System 

C. E. Ekuma, 1 ’ 2 ’E] S.-X. Yang, 1,2 H. Terletska, 1,2 K.-M. 

Tam, 1,2 N. S. Vidhyadhiraja, 3 J. Moreno, 1,2 and M. Jarrell 1, 2, [j] 

1 Department of Physics & Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803, USA 
2 Center for Computation and Technology, Louisiana State University, Baton Rouge, Louisiana 70803, USA 
3 Theoretical Sciences Unit, Jawaharlal Nehru Center for Advanced Scientific Research, Bangalore, 560064, India 

The interplay of interactions and disorder is studied using the Anderson-Hubbard model within 
the typical medium dynamical cluster approximation. Treating the interacting, non-local cluster 
self-energy (E c [G]{i,j ^ *)) up to second order in the perturbation expansion of interactions, U 2 , 
with a systematic incorporation of non-local spatial correlations and diagonal disorder, we explore 
the initial effects of electron interactions (U) in three dimensions. We find that the critical disorder 
strength {WY), required to localize all states, increases with increasing U ; implying that the metallic 
phase is stabilized by interactions. Using our results, we predict a soft pseudogap at the intermediate 
W close to W,Y and demonstrate that the mobility edge (w c ) is preserved as long as the chemical 
potential, /x, is at or beyond the mobility edge energy. 
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Introduction The metal-insulator transition (MIT) 
driven by random impurity has been an important topic 
in physics since the pioneer work by Anderson llj] . A sig¬ 
nificant advance in the MIT theory is achieved by study¬ 
ing it in the context of critical phenomena. Concepts 
from scaling, renormalization group (RG), and random 
matrix theory are used to understand the mechanism of 
localization at different dimensions for different symme¬ 
try classes n®. It has been demonstrated that an in¬ 
finitesimal amount of disorder can lead to localization for 
the models in the orthogonal class at lower (one and two) 
dimensions, whereas there is a MIT for three dimensions 
(3D) 0. In 3D, a sharp mobility edge separating local¬ 
ized and delocalized states develop as disorder strength 
increases @. 


While the MIT of non-interacting systems by now is 
fairly well understood HSU, earlier studies suggested 
that interaction could play an important role in the 
MIT Q. Over the last few decades, experimental works 
ranging from doped semiconductors Si US, perovskite 
compounds jutl - flil ] ~) . to cold atoms in optical lattices fl6l— 
0 have highlighted the importance of the interplay of 
disorder {W) [J 0 0, 0 and interactions (U) 0. 

At the Fermi level, Altshuler-Aronov (0 showed that 
interactions can induce a square-root and logarithmic sin¬ 
gularity in two and three dimensions, respectively, while 
Efros-Shklovskii demonstrated the Coulomb gap [211. 
Field theory perturbative RG method and diagrammatic 
theory which go beyond the Hartree-Fock approxima¬ 
tions have suggested a metallic state for two dimen¬ 
sions [ 22 ], [0. The recent RG work by Finkelstein 
and co-workers has further indicated the possibility of 
a MIT for a model with degenerate valleys 24[, the va¬ 
lidity of which was confirmed through experiments in Si- 
MOSFETs 0[2i]. 

In this letter, we focus on the system with weak local 
interactions on disorder systems in 3D. Our approach is 


an extension of the recently developed typical medium 
dynamical cluster approximation (TMDCA), which has 
shown to be highly successful in describing the Ander¬ 
son localization transition (ALT) for the non-interacting 
systems 0|. The typical medium approaches assume 
that the typical density of states (TDoS), when appro¬ 
priately defined, acts as the “proper” order parameter 
for the ALT. Such an assumption is well justified not 
only for the non-interacting case l2(il [2sl but also in the 
presence of interactions, as shown experimentally0, [0. 
Th e ty pical medium theory (TMT) of Dobrosavljevic et 
al [27| is a special case of the TMDCA when the clus¬ 
ter size N c = 1. Even though the TMT cannot include 
weak localization effects due to coherent backscattering, 
it still does qualitatively predict a disorder-driven ALT, 
and hence incorporates ’strong localization’ effects. The 
TMDCA incorporates non-local effects via systematic fi¬ 
nite cluster increment and achieves almost perfect agree¬ 
ment with numerical exact calculations. The extension 
of the TMT to finite interactions show that interactions 
screen the disorder 3032J. In this letter, we show that 
such a conclusion is robust in the thermodynamic limit 
through increasing cluster size calculations. 


While there have been significant efforts to understand 
the combined effect of disorder and interactions on the 
local density of states close to the Fermi level, the band 
edges have received scant attention. Specifically, the ef¬ 
fect of weak interactions on the mobility edge has not 
been discussed thus far. We are particularly interested 
in the evolution of the mobility edge under the influ¬ 
ence of the Hubbard interaction for spin—1/2 system. 
The transition between metal-the Fermi liquid phase; 
and insulator-the Anderson localized phase is discussed, 
whereas the possibility of the Mott insulator is excluded 
in this study, as only short range, weak interactions will 
be considered. 


The main result of this letter is that for /1 < w e , 
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arbitrary small interactions lead to the masking of the 
sharp mobility edge that separates localized and ex¬ 
tended states in the non-interacting regime below the 
critical disorder strength W^ =0 . Thus, interactions can 
radically modify the spectrum of a non-interacting sys¬ 
tem at the band edges, i.e., in the ‘localized band’. How¬ 
ever, when the chemical potential (/i) is at or above the 
mobility edge energy (i.e., p ^ w e ), the well-defined lo¬ 
calization edge is restored. Nevertheless, unlike the non¬ 
interacting systems where the TDoS just shifts rigidly 
as one scans through p, in the presence of interactions, 
there is a non-trivial decrease of the TDoS vis-a-vis the 
change in the filling. Further, we find a soft-pseudogap 
at intermediate W just below WjS. 

Method The Anderson-Hubbard model (AHM) is 
a model for studying the interplay between electron- 
electron interactions and disorder. The Hamiltonian for 
this model is 

H = - Uj(cl a c ja +hM.)+^(Vi-v)ni a +Uj2n it n zi . 

(ij)cr icr i 

(i) 

The first term describes the hopping of elect rons on the 
lattice, cj^) is the creation (annihilation) operator of 
an electron on site i with spin <r, rii = c\ q is the number 
operator, tij = t is the hopping matrix element between 
nearest-neighbor sites. The second term represents the 
disorder part which is modeled by a local potential V) 
randomly distributed according to a probability distri¬ 
bution P(Vi), p is the chemical potential. The last term 
describes the Coulomb repulsion between two electrons 
occupying site i. We set 4t = 1 as the energy unit and 
use a “box” distribution with P(V)) = 2 W@(W — |V)|), 
where 0(x) is the Heaviside step function. We use the 
short-hand notation: (...) = J dViP{Vi){...) for disorder 
averaging. 

Our focus is on the single-particle Green function and 
the associated density of states. To obtain these for 
the AHM ©, we modify the TMDCA to treat both 
disorder and interactions. Here, an initial guess for 
the hybridization function (T(K,w) = 310 -2 ) is used 
to form the cluster-excluded Green function Q(K,oj) = 
(w — P(K,w) — ck + p) 1 , where £k is the coarse¬ 
grained bare dispersion. Q( K,w) is then Fourier trans¬ 
formed to form the real space Green function, t/„ jm = 
fy(K) exp(*K-(R„ — R m )) and then for a given disor¬ 
der configuration V, we may calculate the cluster Green 
function G C (V) = (fy _1 — V’) -1 . 

Utilizing G C (V), we then calculate the Hartree- 
corrected cluster Green function C/“ 1 (U, U) =G' C (U)“ 1 + 
£d(U) (where £d(U) = p — Uiii/2 and rii = 

— 1 /7r f® oo ^G c (i,i,u))du! is the site occupancy at zero 
temperature, T = 0). Both Q and rii are converged and 
then used to compute the second-order diagram shown 
in Fig. [I] We note that rii obtained at C? c -level is nu¬ 
merically the same as using the full Green function since 


rii is self-consistent at the TMDCA level. This also en¬ 
ables the incorporation of crossing diagrams (for N c > 1) 
from both disorder and interactions at equal footing and 
it is computationally cheaper (for N c = 1, it is ~ 8 times 
cheaper), enabling simulation of large systems. 

Here, we choose the chemical potential p = p + U /2 
to enable simulations both at and away from half-filling. 
Thus, the full self-energy due to interactions is then 
w) = Vi^[Q] + Ec SOPT ^ [f?], where the first term 
is the static Hartree correction and the second term is 
the non-local second-order perturbation theory (SOPT) 
contribution. We note that the computational cost grows 
exponentially with each order of the perturbation series 
making it numerically prohibitive to include more dia¬ 
grams. However, since our focus is on the weak inter¬ 
action regime U/At <C 1, we expect that higher order 
diagrams are suppressed by at least ~ U 3 . 

We have carried out extensive benchmarking of the 
TMDCA-SOPT cluster solver against numerically exact 
quantum Monte-Carlo calculations (UE3 within the dy¬ 
namical cluster approximation (DCA) framework. For 
weak interactions and essentially all disorder strengths, 
the corrections due to perturbation orders higher than 
the second are found to be negligible (for details, see 
Supplemental Material (SM) j42|). 



FIG. 1. The first and 
second-order diagrams 
of the interacting self¬ 
energy between sites i 
and j. 


For a given interaction 
strength U and randomly 
chosen disorder configuration 
V, we calculate the fully 
dressed cluster Green func¬ 
tion g c (v,u) = (g - 1 -V- 
T, Int (U) + U/2)- 1 . With 
G C (K, ui, V, U), we calculate 
the typical density of states 
as 


.V) 

u,V) 

( 2 ) 

following the prescriptions of Ref. |26| , which avoids self¬ 
averaging. The disorder and interaction averaged typi¬ 
cal cluster Green function is obtained using the Hilbert 
transform G£ yp (K,w) = / du>'p1 yp (K, o/)/(w - u/). We 
close the self-consistency loop by calculating the coarse¬ 
grained cluster Green function of the lattice 


P c t yp { K,w) = exp 



(In Pi(w, V)) 


P c ( K,w, 


AC Ei Pit 


V(K u)= [ _ NS{K,e)de _ 

1 J J (Gfj /P (K, w )) -1 + T(K, w) — e + e(K) ’ 

( 3 ) 

where 7 Vq(K, e) is the bare partial density of states. 

Results and Discussion- We start the analysis of our 
results by comparing the algebraic (or average) density 
of states (ADoS) (obtained from the DCA, where the 
algebraic averaging is utilized in the self-consistency) and 
the typical density of states (TDoS) (obtained from the 
TMDCA-SOPT, where the self-consistency environment 
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is defined by a typical medium) for a finite cluster N c = 
38 at various disorder strengths for U = 0.0 and 0.1 at 
half-filling (Figs.EJa) and (b)). 

At weak disorder, W ~ 0.5, the TDoS resembles the 
ADoS. However, for larger W, comparing the U = 0.0 
results (Fig. Ela)) with those of U = 0.1 (Fig. Etb)), a 
noticeable renormalization of the spectrum is observed. 
There is a gradual suppression of the TDoS as the dis¬ 
order strength is increased for both U = 0.0 and 0.1. 
The TDoS at w = 0 is noticeably larger when the U 
is finite. This indicates a delocalizing effect of inter¬ 
actions which is consistent with a real space renormal¬ 
ization group study (43j and has been interpreted as a 
screening of the disorder |30j, 441. For a given disorder 


strength, the band edges at half-filling for the interact¬ 
ing case appear to be identical to that of the U = 0 
spectrum. This seems to imply that the mobility edge is 
preserved when U is turned on. However, this is not the 
case, and this becomes clear upon examining the tails of 
the density of states. 

To explore the effect of weak interactions on the lo¬ 
calization edge of a disordered electron system, we show 
in Fig. [3] the evolution of the TDoS with S = W/W^ 
for various values of U on a linear log plot at various /i. 
Clearly for U = 0, a sharp, well-defined mobility edge is 
observed (see also Fig[2|a)). However, even for a very 
small U = 0.1 (1/30 of the bandwidth), and for both the 
TMT and TMDCA-SOPT, the sharp localization edge is 
replaced by an exponential tail, when /i < ui e . Hence, the 
incorporation of Coulomb interactions in the presence of 
disorder for /i < co e leads to long band tails that are ex¬ 
ponentially decaying. This fingerprint can be understood 
from a Fermi liquid perspective. 

If we inject an electron into a Fermi liquid with an en¬ 
ergy w above the Fermi energy, then, we expect the parti¬ 
cle to experience an inelastic scattering, due to U which is 


W=0.5 


W=1.8 


W=2.0 



FIG. 2. (Color online). Evolution of the ADoS and TDoS at 
various W at U — 0 (a) and 17 = 0.1 (b) for the TMDCA- 
SOPT with N c = 38 for the half-filled Anderson-Hubbard 
model (AHM). 



co <n co 


FIG. 3. (Color online). The evolution of the TDoS of the 
AHM for increasing 17-values for the TMT (N c = 1, (a)) 
and finite clusters (N c = 12 (b) and 38 (c)) at fixed 5 = 
W/Wc = 0.86 on a log-linear plot for increasing /x-values. 
For U = 0.0, we show the plot for fi = 0 only, since changing 
fi only involves a rigid shift of the TDoS. For U > 0, notice the 
systematic disappearance of the exponential tails (indicated 
by arrow) and the non-trivial decrease of the TDoS for the 
finite U (unlike the rigid shift in U = 0) as one approaches 
the mobility edge energy. 


proportional to lo 2 . One factor of ui is due to energy con¬ 
servation and the other to momentum conservation with 
both constrained by the Pauli principle. I.e., the inelas¬ 
tic scattering vanishes as lo — > 0. However, if we apply 
the same logic to an interacting disordered system, then, 
we might expect the edge of the TDoS to be smeared out 
by these inelastic scattering processes, whenever the edge 
energy is above the Fermi energy, but become sharp as 
the edge, approaches it. Though, some argue that this 
reasoning fails for a disordered system, especially for a 
strongly disordered system since a well-defined quasipar¬ 
ticle no longer exists 0 . As a consequence, the concept 
of a mobility edge would not hold and the TDoS should 
have pronounced exponential “tails” even when the Fermi 
energy approaches the top or bottom of the TDoS bands. 
As it is evident from Fig. [3] the sharp mobility edge is 
restored as the mobility edge energy is approached in 
tandem with the Fermi liquid description. 

The smearing of the TDoS edge can further be in¬ 
ferred from the convolutions found in the second order 
(and higher) diagrams (cf. Fig. [lj, which will mix states 
above and below the non-interacting localization edge. 
Consider two such states: one localized and the other 
extended, which are now degenerate due to this mixing. 
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FIG. 4. (Color online), (a) The evolution of the TDoS (at 
w = 0) as a function of the disorder strength W for various 
interactions for N c = 1, 12, and 38 at half-filling. The integral 
J ^sT(K, oj) dK duj vanishes at the same W c as the TDoS for 
a given U (not shown), signifying that the absence of the 
hybridization paths leads to the vanishing of the TDoS. As 
indicated by the arrow, increasing U pushes W c to larger 
values, (b) The interaction dependence of the critical disorder 
W¥ for different cluster sizes N c = 1, 12, and 38 of the AHM 
at half-filling. The unit is fixed by setting 4t = 1. The plot is 
generally in agreement with the results of Ref. [4f|. (c) The 
Wc vs 1/N C on a semi-log plot at U = 0.0 and U = 0.2 for the 
half-filled AHM. Note the systematic and fast convergence of 
Wc with cluster size for both cases. 


- 8=0.80 - 8=0.96 - U=0.20 - U=0.30 

- 8=0.88 - 8=0.98 - U=0.25 - U=0.35 



FIG. 5. (Color online). The TDoS vs energy (w) for N c = 12 
( U = 0.0) (a) and N c = 12 {U = 0.2) (b) at various 5 — 
W/W? showing the formation of a pseudogap at intermediate 
W just before H / c c/=0 ' 2 , which is absent when U = 0. (c) 
Shows the TDoS vs oj at a fixed W (close to W^~ 0 ' 2 = 2.50) 
for various U . Note, the data has been scaled with U. (d) 
Same data as in Fig. [5])c) showing the linear dependence of 
the pseudogap on oj. 


Since these states hybridize with each other, both states 
will become extended 

Next, we explore the effect of interactions on the half- 
filled, disorder-driven localization transition. We show 
in Fig. [4] the evolution of the TDoS at the band cen¬ 
ter, oj = 0, for various cluster sizes. The integrated 
escape rate (/ QT(K,oj)dKdoj) (not shown) character¬ 


izes the rate of diffusion of electrons between the impu¬ 
rity/cluster and the typical medium. The vanishing of 
the hybridization paths leads to a localization transition. 
The TDoS vanishes at the same value of Wjf as the in¬ 
tegrated escape rate. 

Figure 0K a) shows that an increase in U from 0.1 to 
0.5 leads to a concomitant increase in . One can say 
loosely that, the zero-temperature effect of correlations 
is an effective reduction in the disorder strength 0i, 
leading to the increase in W c as indicated by the arrow. 
For the TMT (N c = 1), the increases as 1.83, 1.96, 
2.06, 2.16, and 2.25 for U = 0.1 0.5, while for the 

TMDCA (N c = 12), W% increases as 2.34, 2.48, 2.52, 
2.57, and 2.60, and for for the TMDCA ( N c = 38) as 
2.48, 2.59, 2.65, 2.71, and 2.76 for U = 0.1 - 0.5. 


In Figure |U)b), we show the interaction U dependence 
of the critical disorder strength Wjf for N c = 1, 12, and 
38 for the half-filled AHM. For each of the N c , we obtain 
a correlated metal below the lines, and above we have the 
gapless Anderson-Mott insulator. The trend in both the 
single site and finite cluster are alike (i.e., increases 
with increasing U) except for the difference in . The 
almost linear trend observed for the low U is in agreement 
with previous studies (3l|, 0]. Figure HD c) depicts the 
W T ( / as a function of l/N c at U = 0.0 and U = 0.2 
for the half-filled AHM. Note the systematic and fast 
convergence of W c with N c for both cases. 

We further show in Fig.[5]the evolution of the TDoS(w) 
for N c = 12 at C = 0.0 (Fig. Eta)) and 0.2 (Fig. Etb)) 
for various 5 = W/W%. For finite U a soft-pseudogap, 
which is linear in oj (cf. FigEtd)) develops at the Fermi 
energy (note, this is true irrespective of electron filling) at 
intermediate disorder strengths immediately before the 
system becomes localized. In Fig. Etc), we show that the 
pseudogap is robust as a function of U < 1. Noting that 
we have only short-range interaction, this soft-pseudogap 
cannot be attributed to excitonic effects (which are neg¬ 
ligible here) as in the Efros-Shklovskii theory [Uj]. Also, 
since it occurs only in the TDoS and even for N c = 1, it 
cannot be due to the multivalley structure of the energy 
landscape [ 47 J since a single-site cannot generate a mul¬ 
tivalley energy landscape to sustain sets of metastable 
states, and it should be contrasted from the Altshuler- 
Aronov zero-bias anomaly, which is due to weak non¬ 
local interactions and weak disorder j48j. We ascribed 
this soft pseudogap to the same scenario, which causes 
well-defined mobility edge to only exist when /r ^ oj e . 
U suppresses localization and increases the TDoS. How¬ 
ever, near the Fermi energy the phase space for scattering 
by U is drastically reduced leading to the opening of a 
soft pseudogap. Put differently, the pseudogap is due to 
the suppression of inelastic scattering by U due to the 
Pauli principle and energy conservation. It is linear (cf. 
FigEId)), rather than quadratic in w, due to the lack of 
momentum conservation. 

Conclusions- Based on experiment, theory, and simu- 
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lations, there is a growing consensus that the local den¬ 
sity of states in a disordered system develo ps a highly 
skewed |4^|, log-normal distribution 9. [ 29 I |50| with a 
typical value given by the geometric mean that vanishes 
at the localization transition, and hence, acts as an order 
parameter for the ALT. New mean field theories for local¬ 
ization, including the TMT and its cluster extension, the 
TMDCA, have been proposed. In this letter, we extend 
the TMDCA to weakly interacting systems using second 
order perturbation theory. We find that weak local inter¬ 
actions lead to an increase in W c , with the localization 
edge preserved when the chemical potential is at or above 
the mobility edge energy. For finite U we observe a soft- 
pseudogap for values of the disorder strength just above 
W?. 
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DETAILS OF THE FORMALISM 

In our formulation of both the dynamical clus¬ 
ter approximation (DCA) and typical medium DCA 
(TMDCA), we treat both disorder and interactions using 
the Anderson-Hubbard Hamiltonian 

H = ~ X t ij(4* C j*+ h - C -)+^( V i~P) rl icr + U J2 rli I rli I- 

(ij)cr ia i 

(!) 

where symbols have their usual meanings are described 
in the main text. For the TMDCA, instead of using 
the conventional Matsubara frequency approach, we re¬ 
formulate our formalism in real frequency. This facili¬ 
tates the analysis of the zero-temperature physics, where 


Matsubara-frequency-based approaches may not be ade¬ 
quate. It also avoids the difficulty of analytical continu¬ 
ation of the major observables, the Green function, and 
the self-energy, from the Matsubara-frequency to real fre¬ 
quency. To achieve this, we recall that the second order 
perturbation theory self-energy is the expansion up to the 
second order in U ( 0[U 2 ]) around the Hartree-corrected 
host propagator given as 

S^(i,j» = Ef[e] + E( SOPT )[(?] (2) 

where i and j are site indices. The first term = 

Urii/2 is the Hartree term, while the second term is the 
density-density term which in the Matsubara frequency 
is 


S (50PT)^ 


lim 

iuj— >■ 0 ;+ 


u 2 

TVf/3 2 


y, G( K + Q, iuj n + iv m )G( P + Q ,iui p + iv m )Q(P,iu p ) 

m,p, P.Q 


( 3 ) 


where /? is the inverse temperature. We note further that 
the analytic continuation process may miss important 
features of the observable being studied if not done care¬ 
fully aside its inability to study zero-temperature physics, 
which coincidentally, is the regime we are interested in 
the present study. Thus, to avoid such analytic continua¬ 
tion, and since our cluster problem is solved in real space, 


more also, for the fact that it is numerically more advan¬ 
tageous to work in the real frequency than in the Matsub¬ 
ara frequency, we convert the Matsubara sums (Eq. [3j to 
real frequency integrals using the spectral representation: 
(?(«*>) = / dep(G(e))/(iu) — e), where p{G(e)) = -±SsG(e). 
Defining p^(i,j,u) = the second 

order term of the self-energy in real frequency is 


Pe(w) = U 2 J de 1 de 2 pg(ei)pg(uj - ei + e 2 )pg(e 2) x [n/(-ei)n/(e 2 )n/(-w + ei - e 2 ) + n/(ei)n/(-e 2 )n/(w - ei + e 2 )] 

( 4 ) 


where site labels have been suppressed, pg = —1/n 3(7, 
and rif = l/(e /3e + 1) is the Fermi function. We note 
that pg vanishes only for |w| > 3H/2 [l|, where B is 
the full bandwidth (12i = 3 in our unit). The real part 
of the second order term of the interacting self-energy, 
T,ff >PT on each cluster site is obtained via the 


Hilbert transform 


xi 0p y,j, w ) = 


dui 


,px(i,j,u') 


( 5 ) 


Since l/(x + *0 + ) = V(l/x) — iit5(x) (where ‘V' is the 
principal value of the integral), the non-local SOPT self- 
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energy is then 

£f OPT (y, j,u) = T,^ OPT (i,j,uj) - inp^(i,j, w). (6) 

We note that Eq. [4] scales as an O [iV 3 ], where N is the 
number of the grid points used for the integration. This 
N 3 process is dramatically reduced to scale logarithmi¬ 
cally as N In N using fast Fourier transformation 0] • 


BENCHMARKING THE TMDCA-SOPT 

To validate our developed method and benchmark its 
suitability for studying the Anderson-Hubbard model, we 
compare our results for the Anderson-Hubbard model at 
half-filling with results from the dynamical cluster ap¬ 
proximation (DCA) using the continuous-time quantum 
Monte Carlo method (CTQMC) as the cluster 

solver. The quantum Monte Carlo (QMC) methods are 
powerful tools that enable controlled calculations of the 
properties of large quantum many-particle systems. De¬ 
tails of the CTQMC formalisms are well described in the 
literature (see for e.g., Refs. 0-i) as such, and we will not 
attempt to give a detailed description of the algorithms 
here but will just give an overview which enables us to 
compare the QMC results with our Hilbert transformed 
imaginary frequency data. In the CTQMC, to avoid any 
possible time discretization error in the imaginary time 
axis, we adopt the recent improvements in the QMC al¬ 
gorithms in the continuous imaginary time [3, 0, 0, 0. 
This significantly improves the quality of the data. 

To make the comparison, we convert our real frequency 
data to Matsubara frequency using the Hilbert transfor¬ 
mation, and obtain the local Green function and self¬ 
energy, respectively, as 


Gloc{i^ri) 


i r ^3^) 

7T J_ 00 iuJ n -UJ ’ 


£zoc(^n) 


i r 

7T J_ 00 iuj n - U 


(7a) 


(7b) 


where = (2 n + 1)7 t//3 with /3 = 1/ksT and n £ Z is 
the fermionic Matsubara frequency, and ks (=1) is the 
universal Boltzmann’s constant. 


Limit of Zero Disorder: Hubbard Model 

As a natural consequence, the self-energy of the Hub¬ 
bard model of the numerically exact DCA-CTQMC con¬ 
tains the most vital information for benchmarking with 
a mean-field theory like the TMDCA-SOPT. We show 
in Figs. [T]( a) and [Tfb) the comparison plots of our 
TMDCA-SOPT (using Eq. I7bl) as compared with the 



FIG. 1. (Color online). The comparison of the single¬ 
particle quantities from TMDCA-SOPT and DCA-CTQMC 
at T = 0.025 for the half-filled Hubbard model (zero disorder 
and the finite Coulomb interaction U > 0). The compari¬ 
son of the imaginary part of the local self-energy for N c = 1 
(a) and N c = 14 (b), respectively, on a log-log scale. Also 
shown is the corresponding comparison of the local Green 
function for N c = 1 (c) and N c = 14 (d), respectively on 
a linear scale. Smaller values (U < 0.3 for N c = 1) and 
(U < 0.4 for N c = 14) are not shown on the local Green 
functions plot for easy readability as they are too close to 
each other. In both cases, the dotted lines with symbols are 
for the TMDCA-SOPT while the solid lines depict the DCA- 
CTQMC results. Further, for the self-energy plots, smaller 
interaction strengths (U < 0.5) are on top of each other and 
as such, may appear indistinguishable. 


DCA-CTQMC data for the N c = 1 and finite cluster, 
14, respectively, at various values of the interaction. As 
it is evident from the plots, our data benchmarks well up 
to U ~ 2.25 for N c = 1 and U ~ 0.75 for N c = 14 show¬ 
ing both quantitative and qualitative agreement between 
our TMDCA-SOPT data and the DCA-CTQMC results. 
The remarkable exact agreement between the TMDCA- 
SOPT and DCA-CTQMC for low U (< 0.75), at all fre¬ 
quencies for both single-site and finite cluster shows that 
our self-energy has the correct behavior and as such, en¬ 
sures that the conclusions arrived in the main paper are 
numerically correct. At least for the smaller U-values 
(which is the regime we are interested in), the good agree¬ 
ment with DCA-CTQMC further shows that the pertur¬ 
bation expansion of the self-energy up to 0[(f//H) 2 ] in 
U should capture all the dominant quantum fluctuations 
guaranteeing that quantitatively correct results are ob¬ 
tained. Hence, the results we presented in the main pa¬ 
per which are for small interacting disordered electron 
system are accurate to within the computational accu¬ 
racy of our formalism. 

We further show in Figs. □Jc) and □Id), the plot of 
the imaginary local Green function of our data (trans- 
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FIG. 2. (Color online). The comparison of the self-energy from TMDCA-SOPT and DCA-CTQMC at T = 0.025 for the 
half-filled Anderson-Hubbard model. From left to right, the Coulomb interaction U is fixed at 0.1, 0.2, and 0.5, respectively. 
Here, the solid lines are for the DCA-CTQMC while the dash lines with closed symbols are for the TMDCA-SOPT. 


formed from real to Matsubara frequency using Eq. I7al) 
as compared to the DCA-CTQMC data for T = 0.025 
and W = 0 for the N c = 1 and N c = 14, respectively. 
Again, as a confirmation of the good agreement in the 
local self-energy, the local Green function from the two 
methods is numerically the same especially for the small 
[/-values. They are practically on top of each other up 
to U « 1.50 for both N c = 1 and 14. 


Finite Disorder and Interaction: Anderson-Hubbard 
Model 

To further check the applicability of the expansion of 
the interacting self-energy in powers of U up to the sec¬ 
ond order for the study of interacting disordered electron 
system, we further benchmark our developed method for 
the Anderson-Hubbard model with the DCA-CTQMC. 
This again becomes imperative as we are not aware of a 
prior benchmarking especially for the finite cluster. 

We show in Figs. [2] the comparison of our developed 
method for fixed interaction strengths (U = 0.1, 0.2, and 
0.5, respectively) at various disorder strengths, respec¬ 
tively (using Ea l7bl) . with the imaginary frequency data 
of DCA-CTQMC. As it is evident from the plots, there 
is almost a perfect agreement between our Hilbert trans¬ 
formed data with the DCA-CTQMC results even for rel¬ 
atively large value of the interaction strengths (U ~ 0.5) 
for all the disorder strengths up to the localization tran¬ 
sitions. This further confirms that our truncation of 
the perturbation series expansion of the interacting self¬ 
energy at (9[(f//.B) 2 ] (at least in the weak interaction 
regime) is enough to account for the quantum fluctu¬ 
ations in the typical environment. One can thus af¬ 
firm that at least, within the weak interaction regime 
(U/At -C 1) in a disordered electron system that higher 
order diagrams in the expansion of the E* nt are sup¬ 
pressed by at least ~ ( U 3 ). 


EXPLORING THE MOBILITY EDGE 

We show in Fig. [3l the plot of the typical density of 
states on a log-linear scale at half-filling for the TMT 
(N c = 1) and for the finite clusters (N c = 12 and 38) for 
increasing {/-values at different disorder strengths for 3D. 
Clearly for U = 0, the slopes of the TDoS for both the 


N=1 N =12 N =38 



co <n co 


FIG. 3. (Color online). The typical density of states of the 
Anderson-Hubbard model at half-filling for increasing U-value 
for the TMT (N c = 1) and finite clusters (N c = 12 and 38) 
at various disorder strengths on a log-linear plot. Observe 
that the U = 0 results have infinite slopes signifying the lo¬ 
calization edge while for all the interactions, even for a very 
small U = 0.1 which is l/30th the bandwidth, the slopes re¬ 
main finite out to the band edge of the ADoS (see the main 
text) consistent with the mixing of the localized and extended 
states with different energies leading to the suppression of the 
mobility edge (see main text). 
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TMT (N c = 1) and TMDCA (N c = 12 and 38) becomes 
infinity signifying the existence of a sharp, well-defined 
mobility edge. However, for a very small U = 0.1 (1/30 
of the bandwidth), and for both the TMT and TMDCA- 
SOPT, the slope is evidently finite up to high frequencies 
signifying that the sharp well-defined localization edge is 
replaced by an exponentially fast cross-over. Hence, the 
incorporation of Coulomb interactions in the presence of 
disorder at half-filling leads to long band tails that are 
exponentially decaying. 

One can argue that since the tails are far away from the 
chemical potential, /r when /r = 0, they do not generally 
participate in transport since only states in the proxim¬ 
ity of the Fermi energy is excited. This can be attributed 
to the fact that at half-filling with finite U, interaction 
induces the mixing of the localized and extended states 
with different energies of the states at the mobility edge 
leading to the suppression of the sharp, well-defined mo¬ 
bility edge boundary with the emergence of long tails 
with exponentially fast crossover. This is consistent with 
the delocalization nature of U at ui = 0 states at half¬ 
filling leading to the increase in W c . As demonstrated 
in the main text, the mobility edge is well defined even 
in the presence of interactions as long as the chemical 
potential is located at or beyond the mobility edge en¬ 
ergy of U = 0. Hence, Fig. [3] essentially shows that at 
half-filling, due to the mixing of states induced by the 
interaction, the mobility edge energies develop tails with 
exponentially fast crossover which masks its detection. 
However, if we move the chemical potential energy at or 
outside the mobility edge energy (see the main text), the 
sharp, well-defined mobility edge is restored. 

The masking of the localization edge in the TDoS due 
to the mixing of the the states induced by the small in¬ 
teraction can further be confirmed from the convolutions 


point of view. This can be seen from the convolutions 
found in the second order (and higher) diagrams. These 
convolutions will mix the states above and below the non¬ 
interacting localization edge. Consider two such states: 
one localized and the other extended, which are now 
degenerate due to this mixing. Since the elastic disor¬ 
der scattering causes these states to hybridize with each 
other, both states will become extended |lf)j]- 
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